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We present the results of numerical analysis of a model of normal zone propagation in coated 
conductors. The main emphasis is on the effects of increased contact resistance between the super- 
conducting film and the stabilizer on the speed of normal zone propagation, the maximum temper- 
ature rise inside the normal zone, and the stability margins. We show that with increasing contact 
resistance the speed of normal zone propagation increases, the maximum temperature inside the 
normal zone decreases, and stability margins shrink. This may have an overall beneficial effect on 
quench protection quality of coated conductors. We also briefiy discuss the propagation of solitons 
and development of the temperature modulation along the wire. 
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I. INTRODUCTION 

Quench protection of large scale devices, such as mag- 
nets and cables, based on coated conductors has emerged 
as one of the major unresolved obstacles in their imple- 
mentation. High operating temperature and, correspond- 
ingly, relatively large heat capacity make coated conduc- 
tors very stable in comparison to the conventional low 
temperature superconductors. However, the side effect of 
this positive quality is that when a normal zone does nu- 
cleate it expands very slowly. The potential drop across 
a short normal section of a long conductor is difficult to 
detect and in adiabatic or nearly adiabatic conditions the 
temperature of this section may rise above the safe limit 
resulting in irreversible damage to the whole coil or cable 
strand. 

This article presents the results of a numerical analysis 
of a model of normal zone propagation (NZP) specialized 
to the architecture of the state-of-the-art coated conduc- 
tors. Its main purpose is to elucidate the effects of the 
interfacial resistance (contact resistance) between the su- 
perconducting YBa2Cu307_2: (YBCO) film and copper 
stabilizer on stability and speed of NZP. The interest to 
this problem arose initially from an effort to understand 
some peculiar effects that accompany quench in coated 
conductors [TJ[2]. It seems clear now that these phenom- 
ena result from a large resistance between the YBCO 
film and a metal substrate |3j |4j. This understanding 
has lead to realization that increasing the contact resis- 
tance between the YBCO film and copper stabilizer may 
have beneficial effect on the speed of the normal zone 
propagation [5l[6l[7]. 

The effects of a large contact resistance between the 
stabilizer and conventional low T^. superconductors have 
been studied extensively in the past [S] [HI [111 IHl El]- 
However, the idea of tailoring the properties of the su- 
perconducting wires by increasing the contact resistance 
has not been adopted to wider use. This option of con- 



ductor design had lain dormant for many years — a solu- 
tion in wait of a problem. Perhaps, coated conductors 
present just such a problem. Increasing the contact re- 
sistance does make the conductor less stable. However, 
since coated conductors are inherently much more stable 
than the low superconductors, the reduction of the 
stability margins accompanied by increasing the speed of 
NZP may allow to develop coated conductors overall bet- 
ter suited for large scale applications than their current 
version with a minimized contact resistance. 

Here we will discuss the NZP in a straight coated con- 
ductor cooled from the surface. This model more closely 
describes the typical conditions in the experiments such 
as [TJ 12] or in a superconducting cable, rather than in 
a pancake coil. Correspondingly, the concrete example 
presented below is based on the operating temperature 
equal to 65 K and the values of the material parameters 
in the temperature range 65 to 77 K. Our main conclu- 
sions that the increased contact resistance leads to in- 
creased NZP speed and reduced stability margins will 
remain qualitatively valid for any type of application. 
However, in order to adequately describe a pancake coil 
one needs to take into account the heat transfer in the 
radial direction, between the turns, in addition to the 
lateral heat flux which we consider here. The problem of 
NZP in a pancake coil will be addressed elsewhere. 

This paper is structured as follows. In Section II we 
briefly formulate the model of a quench in coated con- 
ductor based on conditions of energy and charge conser- 
vation. The model is similar to that in [SJ [7], but with 
a modified constituent relationship between the electric 
field and current density in the superconductor. How- 
ever, in [5l |7] the numerical solutions were obtained us- 
ing an approximation that is strictly valid only in the 
limit of a small contact resistance (see detailed expla- 
nation in [7]). Therefore, the results of that analysis 
could not be extended to the conditions when the con- 
tact resistance is arbitrarily large. Moreover, legitimate 
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questions remained as to what extent the results of (SJ |7] 
were influenced by the approximation used to solve the 
problem, rather than the physics of the phenomenon. 
The solutions of the system of coupled nonlinear partial 
differential equations were obtained rigorously without 
approximations. Qualitatively, the current results are 
similar to those in [5;. The propagation speed increases 
and stability margins decrease with increasing contact 
resistance. Quantitatively however, there are substan- 
tial differences in the rate of change of these character- 
istics with the value of contact resistance. The rigor- 
ous solutions also confirm the emergence of the dissipa- 
tive structures — spontaneous modulation of temperature 
along the conductor — reported previously [7j. 



II. MODEL 




FIG. 1: A sketch of the cross-section of coated conductor (not 
to scale). Shown are four constituencies we have taken into 
account in this model. These are copper stabilizer (1), metal 
substrate (2), a thin superconducting film deposited on top of 
the substrate (4), and the interfacial resistive layer (contact 
resistance) that segregates the superconducting film from the 
stabilizer (3). 

Coated conductors [T3J |T3] are manufactured in the 
form of a tape in which the superconducting YBCO film 
of about 1 |J.m thick is deposited on a buffered flexible 
metal substrate (e.g. Ni-W alloy, Hastelloy or stainless 
steel). A copper stabilizer is either soldered or electro- 
plated on top of the YBCO fllm. See Figure [l] The 
standard width of such a tape-like wire is 4 mm, the 
thickness, about evenly divided between the substrate 
and stabilizer, is close to 100 |i.m. 

In [7] it was shown that the 3D equation of heat con- 
duction in a thin tape-like composite wire can be reduced 
to a 2D (planar) or ID (linear) model if the heat flux 
from the surface does not greatly exceed 1 W/cm^. In 
this case the temperature variation across the thickness 
of the tape constitutes a fraction of a degree and is much 
smaller than the variation of temperature along the wire. 
The ID approximation is valid as long as the thermal dif- 
fusion length is greater or comparable to the conductor 
width. In coated conductors this condition is met. The 
ID (in-plane) heat conduction equation for the coated 
conductor has the form [7] 
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d_ 
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K 
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^Q-2KoiT-To). (1) 



Here C = Cidi + C2d2 is the combined heat capacity, 
K = Kidi + K2d2 is the effective thermal conductivity. 
The subscripts 1 and 2 refer to the stabilizer and sub- 
strate, respectively. The thicknesses of the stabilizer and 
substrate are denoted as di and d2- Q = /^^^ q{z)dz is 
the density of the internal heat sources integrated over 
the thickness of the wire. Kq is the heat transfer coeffi- 
cient across the insulation on the surface of the wire and 
To is the ambient temperature. 

The instantaneous redistribution of current between 
the superconducting fllm and stabilizer is determined by 
the condition of charge conservation 



dJi 
dx 



Vi-Vs 
R 



(2) 



where Ji is the linear density of current (A/cm) flowing 
through the stabilizer. Vi and Vs are the local electric 
potentials of the stabilizer and superconductor, respec- 
tively and R [fi cm^] is the resistance of the unit area 
of the interface (contact resistance). This condition can 
also be used in the form 



ox \ ox 
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Here Ei = —dVi/dx and Es — —dVs/dx are the electric 
flelds in the stabilizer and superconductor, respectively, 
and we do not assume that R is uniform. 

The integrated area density of heat sources in ([T|) takes 
the form 
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(4) 



where = J — Ji is the density of current in the super- 
conductor and J (constant) is the total transport current 
density in the coated conductor. 

The constituent relationship for a superconductor can 
be presented in many forms. Here we will use the one 
from [15j 

E,{J,) = i?„ Jo In (1 + exp{( J, - J,)/M) . (5) 

Here Jc{T) is the critical current and i?„ and Jq are phe- 
nomenological parameters. They can be determined from 
the limiting cases. It is customary to deflne the critical 
current by the condition that at Jg — the electric fleld 
in the superconductor is equal to E^ = \ \^^^ /cn\. Thus, 



i?„Jo = ^o/ln2. 



(6) 



On the other hand, when Js substantially exceeds Jc the 
electric fleld is determined by the flux flow 



Es{Js) -J-i-^s — Jc), 



(7) 



where ps is the normal state resistivity and ds the thick- 
ness of the superconducting fllm. Thus, 



R„ 



(8) 
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and 



Jo 



Eods 



(9) 



In coated conductors dg ^ I [un and at 100 K the value 
of ps ~ 100 [iTl cm, so that 



Jo - 10"*^ A/cm. 



(10) 



Thus, for all practical purposes ([5| can be used in the 
piecewise form 



Rn{Js-Jc), iiJs>Jc 

0, if Js < Jc 



(11) 



Hereafter, we will adopt a linear temperature dependence 
of J, inifMl: 



Jc^a{T,~T); T < T,. (12) 
At T > Tc the superconductor has ohmic resistance 



(13) 



For the stabilizer, the conventional Ohmic relationship 
will suffice at all temperatures: 



E,{,h) = ^Ji. 
di 



(14) 



In the normal state the resistance of YBCO film is much 
greater than that of the stabilizer. 



ds di 



(15) 



The final step in formulating this model is to 
present ([T|) and ([3| in the dimensionless form. The cur- 
rent sharing temperature Ti is defined by the condition 
Jc{Ti) — J. Let us introduce a dimensionless tempera- 
ture 6 



e 



T -Ti 



Then, ( 12 1 takes the form 

J, = J(l-0). 



(16) 



(17) 



Let us introduce a fraction of the total current that flows 
through the stabilizer 

Ji ^ Ju; Js = J{l~u); 0<u< 1. (18) 



For 9 < 1 equation (|11|) takes the form 
E, 

For 6* > 1 



RnJ{9-u), u<e 
0, u>6 



^J(l-u); 



(19) 



(20) 



To avoid an unphysical discontinuity at 61 = 1, we will 
consider i?„ = ps/ds 
compact form 



Then, (|3| can be written in a 
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(22) 



(23) 



is the current transfer length which determines the length 
scale of the current exchange between the superconductor 
and stabilizer [4]. Taking into account Q, the first two 
terms in the right hand side of Q take form 



Pi J 2 
di 



du 
dx 



(24) 



Taking into account ([3|, the last term — losses in the 
superconductor — can be written as follows 



j,£;, = j2(i-u) 



d_ 
dx 



R—\ 
dx I 



(25) 



We will express the distances in units of thermal dif- 
fusion length It and time in units of 7^^, where 



li 



{DtHY'^; j = piJydiCAT. 



(26) 



Here Dt = K/C is the effective thermal diffusivity of 
the conductor, AT = Tc — Ti, and the increment 7 deter- 
mines the characteristic time required for the Joule heat 
generated in the stabilizer to warm the conductor by the 
temperature AT. 

In dimensionless variables, (jlj and (21) take the form 



00 
dr 



d_ 
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with T — "ft and ^ 

_ 2KoATdi 
" PiJ' 



= u — T max [0, min(0, 1) — u] . 



x/It- Here 

'o = (To-Ti)/(T,-Ti) <0. 



(27) 



(28) 



(29) 



Notice that (27 1 does not depend on the specific form 



of the constituent relationship between electric field and 
current density in the superconductor. The specifics of 
the constituent relationship enters only in the charge con- 
servation condition given by ([3| and its dimensionless 
versions (21 1 and (28 1. 
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The relative role of the interface resistance is deter- 
mined by the parameter 



A. Speed of normal zone propagation 



R 

Rq ' 



Rn — 
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di 



KAT 



(30) 



The results will not depend on the value of F as long as 
r ^ 1. Hereafter, for the purpose of numerical calcula- 
tions, we take F = 10^. 



III. RESULTS 

As the first step let us consider the interplay between 
the heat source and cooling power for uniform temper- 
ature {d'^e/di'^ = 89/ di = 0). The analytical solu- 
tion u{9) of ( [28| can be used in (27 1 to determine the 
heat source as a function of temperature. The result 
is practically the same (as long as F ^ 1) as that 
shown in Figure 3 in [7j. The system is bistable when 
K < Kc = 1/(1 + l^oD- In this case there are two stable 
uniform modes of operation. One is high temperature 
^max — 9{)-\- and high dissipation and the other with 
low temperature ^min — and zero dissipation. When 
K > Kc (cryostable condition), there is only one stable 
uniform state — the low temperature state with tempera- 
ture 9q. 



The system of equations (27 1 and (28) were solved 



numerically by using an IMEX Crank-Nicolson/Adams- 
Bashforth method in conjunction with a fixed-point 
method to solve the Poisson equation with nonlinear 
source term. The purpose of the numerical solutions 
is to determine the regions of the physical parameters 
{r, K, ^o} that correspond to different types of conductor 
response to initial perturbation, such as the normal zone 
propagation or formation of the dissipative structures. 
We also determine the speed of NZP and the margins of 
stability as the functions of the contact resistance. 

The solutions t) presented below correspond to 
periodic boundary conditions and the initial condition in 
the form of a Gaussian in the center of the conductor 



9{tQ) = {a-9o)c-^"'^'" +9o. 



(31) 



For (28) we used the following boundary conditions: 



u{L) = u{~L) = 1. 



(32) 



This condition means that the current is injected at the 
ends of the conductor into the stabilizer, as is the case 
in real experiments. It is important to emphasize that 
this boundary condition requires that the contact resis- 
tance R was very low at the ends of the conductor. This 
allows the current injected into the stabilizer to transit 
into the superconductor over a very short distance with- 
out generating much heat at the ends. Without a low 
contact resistance at the ends the conductor will become 
unstable when contact resistance exceeds a certain level. 
In numerical calculations we have used R{x) and, corre- 
spondingly, r(^) that are constants everywhere, except at 
the ends of the conductor near ^ = ±_L where r(^) <C 1. 
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FIG. 2: Solutions of ([27| and ([28|) for two different values 
of the contact resistance. Dark (red) color indicates ele- 
vated temperature. The normal zone propagates with con- 
stant speed determined by the slope dx/dt a,t T — Tc {9 ~ 1). 
The time scale here < 7i < 80 and the length of the sample 
— 50 < x/It < 50. The values of the contact resistances de- 
fined by the parameter r^^^ = \/It are indicated. The upper 
scale and the scale on the right show the distance along the 
conductor in centimeters and elapsed time in seconds, respec- 
tively. These values correspond to the specific set of material 
and operating parameters defined as an example in the text 
[It = 0.4 cm and 7 = 6.7 s"^). 

Figure|2]illustrates the growth of the normal zone prop- 
agation, which emerges from the initial Gaussian pro- 
file. The figure illustrates the differences between the 
solutions that correspond to normal zone propagation at 
different values of the contact resistance. In each case 
the cooling constant k = 0.2 and the operating tem- 
perature 9o = —1. According to (17) this corresponds 
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to the conductor operating at 50% capacity, namely 
J° = Jc{To) = 2J. The propagation speed noticeably 
increases with increasing A. It should be noted, that for 
the finite values of the cooling constant k the temper- 
ature behind the propagating front does not depend on 
the contact resistance. In the adiabatic case the situation 
is different and is discussed in the next subsection. 

R cm^) 

128 288 




FIG. 3: The NZP speed in units of Ut as a function of contact 
resistance {X/It = {R/ Ro)^^^) for several values of the cooling 
constant k. For illustration purposes the scale on the right 
and the upper scale show the NZP speed and the values of the 
contact resistance for a set of material parameters described 
in the text {Ut = 3 cm/s, Ro = 8 \iQ cm^). 



The natural unit for NZP speed is 



Ut — ItI = 



diC^AT J 



1/2 



(33) 



In Figure [3] the speed of normal zone propagation ex- 
pressed in units of Ut is shown as a function of contact 
resistance for several values of the cooling constant. For 
X/It > 2 the speed scales approximately with A which 
means that the NZP propagation speed is determined by 
the greater of the two length scales specific to this prob- 
lem: 



V oc max[/T, A]7. 



(34) 



Let us to fiesh out these conclusions using the val- 
ues of the material parameters representative of coated 
conductors [16j . For copper stabilizer we will take 
Ki w (4to5)W/cmK; Ci w l.TJ/cm^K; pi w 
0.2 X 10~^ cm, and di = 40 \mi . For Hastelloy 
substrate we take K2 ~ 7 x 10^^ W/cm K; C'2 ~ 
1.4 J/cm"^ K, and ^2 — 50 ixm. Let us also take the crit- 
ical temperature Tc « 90 K and the operating tempera- 
ture To = 65 K. A reasonable self field value of the crit- 
ical current density Jc(7o) = 300 A/cm, and the trans- 
port current density J = 150 A/cm (the corresponding 



value of 00 = ~1) • Then, the current sharing tempera- 
ture Ti = 77.5 K and AT = - Ti = 12.5 K. 

The effective thermal conductance of the wire is dom- 
inated by the copper stabilizer 

K = Kidi + K2d2 « (1.6 to 2) X 10"^ W/K. (35) 

The combined heat capacity 

C = Cidi + Csda « 1.4 x 10"^ J/cm^ K. (36) 

The corresponding thermal diffusivity and the increment 

Pi J' 



Dt « (1.14 to 1.4) cmVs; 7 = 



diCAT 



6.7 s- 



(37) 



The length scale in this problem is determined by the 
thermal diffusion length (26 1, 



It w (0.4 to 0.45) cm. 



The natural scale of the propagation speed ( 33 1 
Ut ~ 3 cm/s. 



(38) 



(39) 



The contact resistance in currently manufactured coated 
conductors reported in [17j 



R. 



5 X 10"^ n cm^ 



(40) 



is consistent with other data reported in literature. For 
this value of the contact resistance the current exchange 
length A (pSl), 



A w 3 X 10"^ cm < It- 



(41) 



In the limit of low contact resistance, X/It <C 1, the 
propagation speed V is smaller than Ut in agreement 
with experimental findings in which the NZP speed was 
found to be in the range of ( 1 to 2) cm/s [THUT^ PUl . 

The characteristic contact resistance defined by the 



condition X — It, pO]) 
Ro 



di 



^ w 8 X 10"*^ r^cm^ 



(42) 



Thus, in order to achieve a substantial increase in normal 
zone propagation speed, the contact resistance has to be 
increased well over 10 [iH, cm^. 



B. Temperature rise inside normal zone 

The main danger of slowly propagating normal zone 
is that it may remain undetected for an extended period 
of time during which the temperature inside NZ will rise 
above the damage threshold. This is especially true for 
adiabatic or near adiabatic conditions. The voltage drop 
across the normal zone and, correspondingly, the proba- 
bility of its detection increases approximately in propor- 
tion to the length of NZ. Therefore, in order to improve 
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FIG. 4: (a) The maximum temperature 6'niax inside NZ as 
a function of its length ^„ in the adiabatic case (k = 0) for 
three different values of the contact resistance. The scale 
on the right shows the values of T^ax = Ti + ^maxAT with 
Ti = 77.5 K and AT = 12.5 K. The top scale shows the 
length of the NZ in centimeters {It ~ 0.4 cm), (b) The rate 
of temperature increase d6yaa.yi/ d(^n versus contact resistance. 
The scale on the right shows dTma.^/dln = d9mzx/dS^n{^T /It). 
The upper scale shows the values of the contact resistance 
corresponding to i?o = 8 err? in ( 30 1 . 



the quench protection quahty of the conductor the rate 
at which the temperature inside NZ rises with the length 
of NZ has to be made smaller. Obviously, the increas- 
ing speed of NZP accomplishes just that. However, the 
increased contact resistant also increases the amount of 
power dissipation at the front of the propagating NZ. 
This makes it necessary to examine closely how the max- 
imum temperature inside the normal zone changes with 
its length in the worst case scenario of adiabatic NZP 
{n = Q). 

Figure |4]^a) shows how the peak temperature changes 



with the length of NZ Z„ defined as the length of a section 
with T > Tc- In adiabatic conditions the temperature in- 
side NZ increases practically linearly with its length. One 
can clearly see the benefit of increased contact resistance. 
When the length of NZ reaches 40 cm {In/h = 100) the 
peak temperature of the hot spot may reach « 750 K 
in the conductor with low contact resistant X/It ^ 1- 
In the conductors with substantially larger contact resis- 
tance, such that X/It = 4 to 8, the peak temperature is 
substantially lower for the same length of NZ. Of course, 
the calculated values of temperature here are given only 
for the purpose of comparison because we do not take 
into account the changes with temperature in resistivity 
and other material parameters. 

A more detailed picture of the effect that the contact 
resistance has on the temperature inside NZ is given in 
Figure |4|^b). The figure shows the rate dT^^^^/dln at 
which the peak temperature increases with the length of 
NZ (the slope of the curves in Figure [4]j a)) as a function 
of the contact resistance expressed as X/It- This rate al- 
lows us to estimate Tmax for an arbitrary length of the NZ 
at a given value of the contact resistance. For X/It ^ 1 
the temperature of the hot spot increases approximately 
by 18.5 K/cm of the NZ length. In a conductor with 
X/It ~ 8, the rate of the temperature increase is less 
than 5 K/cm. The inset to Figure [4]^b) shows the con- 
secutive temperature profiles of the adiabatic NZP for 
X/It = 4. 



C. Stability margins 

An ability of a coated conductor to dissipate a certain 
amount of heat deposited by an external source with- 
out triggering normal zone propagation generally declines 
with increasing contact resistance. We illustrate the re- 
duction in stability margins by using the initial condition 
in the form of (31 ) with fixed width 5 — \ A and variable 



peak temperature Tp determined by the parameter 



(43) 



Physically, this corresponds to a rapid injection of a cer- 
tain amount of heat into a small section of the conductor. 
For a given value of the contact resistance we determine 
the value of a above which the initial temperature profile 
gives rise to NZP. For smaller values of the peak tem- 
perature the initial profile dissipates without triggering 
NZP. In Figure |5] the stability boundaries are shown for 
several values of the cooling constant. The phase space 
{Tp, X/It} below the respective curve corresponds to the 
range of stability, perturbations above the curve are un- 
stable. For illustration purpose the scale on the right 
shows the peak temperature in absolute units for the set 
of material and operating parameters described above 
(To = 65 K,Ti = 77.5 K,Tc = 90 K). 

It is obvious that the stability margins precipitously 
decline with increasing contact resistance and become 
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FIG. 5: Stability boundaries as functions of X/It for different 
levels of cooling. The scale on the left shows universal dimen- 
sionless units of the peak temperature, Eq. (43). The scale on 
the right shows the peak temperature in absolute units for the 
specific set of material and operational parameters described 
in the text, Tp = Ti + aAT. When the peak temperature of 
the initial temperature profile exceeds the limit shown by the 
curves, the normal zone starts to propagate. 



less dependent on the cooling conditions. However, our 
analysis of this model shows that the stability of the ini- 
tial perturbation depends also on the width of the pertur- 
bation. The greater the width of the initial temperature 
profile, the more stable the conductor is with respect to 
that perturbation. When the width of the initial temper- 
ature profile is greater or comparable to A its stability is 
about the same as that of the conductor with X It- 
Second, even in the worst case scenario the conductor re- 
mains stable as long as the peak temperature is below Ti 
(JciTp) > J) because no current will be diverted into the 
stabilizer and no heat will be generated (at least within 
the constituent relationship (11) we have adopted here). 
For the set of parameters described above the minimum 
amount of heat per unit area that can be dissipated with- 
out triggering NZP can be estimated as 



water and sand dunes on the ground, to a planetary size 
phenomenon like the north pole hexagon — a long-lived 
feature of the atmosphere of Saturn [53] . 

The spontaneous temperature and critical current 
modulations along the superconducting wire take place 
under conditions of strong cooling, above or very near 
cryostability condition k > k^, = 1/(1 + l^'ol))- At this 
level of cooling the maximum amount of heat generated 
in the stabilizer is not large enough to maintain temper- 
ature above the critical. However, if the contact resis- 
tance is large enough, the additional heat generated in 
the interface by the current passing between the super- 
conductor and stabilizer may be sufficient to sustain a 
modulated temperature profile. 

A detailed analysis of different scenarios and types 
of modulation will be published elsewhere. Figure [6] 
presents an example of spatio-temporal development of 
the initial Gaussian perturbation in cryostable regime, 
K = 0.6, 6*0 = -1, and X/h = 5.48 {R/Rq = 30). The 
color denotes temperatures from the operating temper- 
ature (65 K) to the maximum temperature equal to Tc 
(90 K). 

The evolution proceeds as follows: The initial tem- 
perature peak splits into two twin peaks. Each of those 
peaks, in turn, continue the division, but one of the twins 
produced in each cycle does not survive. The surviving 
twin is shifted from the position of its parent and starts 
the same cycle of division. This results in the directed 
drift of the soliton (hot spot) along the wire. Eventu- 
ally, the two moving solitons are separated by a distance 
large enough to even-out the chances of survival of both 
twins and one of the cycles of division ends successfully, 
doubling the number of moving hot spots, etc. The peak 
temperature of the hot spots is close to the critical tem- 
perature. Figure|6jb) shows the temperature distribution 
along the conductor at a moment jt = 150. This pro- 
file corresponds to the cut shown by the dotted line in 
Figure [6]^ a). 

It should be noted that a similar type of temperature 
modulation was discussed in O UHl HH HI] , where they 
were called resistive domains. In |llj an experimental 
observation of such a resistive domain was reported. 



C{Ti - To) - 0.17 J/cm^ 



(44) 



IV. SUMMARY 



D. Solitons and dissipative structures 

Besides the two uniform modes of operation discussed 
previously a superconducting current-carrying wire may 
exhibit a more complex behavior that can be classified as 
the formation of dissipative structures [7]. The spatially 
uniform physical systems driven away from thermal equi- 
librium tend to break the translation symmetry and form 
steady or time- variable macroscopic spatial patterns |22j . 
The examples run the gamut from the table-top demon- 
strations of convection cells, to the sand ripples under 



There is a viable option to improve the quench pro- 
tection quality of coated conductors by increasing the 
contact resistance between the superconducting film and 
stabilizer. This increases the normal zone propagation 
speed but decreases the stability margins with respect to 
localized temperature perturbations. A compromise be- 
tween these two requirements can be found which may 
yield a better overall superconducting wire. We should 
emphasize that the contact resistance here is the resis- 
tance to the current exchange between the superconduc- 
tor and stabilizer, not the resistance to the current fiow 
through the stabilizer itself. It is still desirable to have 
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copper stabilizer of substantial thickness in order to min- 
imize the Joule heat generated in the stabilizer. The 
contact resistance can be introduced by various means. 
For example, a very thin film (perhaps 100 nm thick) 
of highly resistive substance can be deposited on top of 
YBCO and then covered with protective silver layer and 
copper stabilizer. 



The model presented here describes the process of NZP 
in a straight conductor where the heat is transferred only 
along the wire. This condition is similar to that in a 
superconducting cable. For a pancake coil our conclu- 
sions about increasing speed of NZP and reduced sta- 
bility margins as the functions of the contact resistance 
remain qualitatively valid. However, a thorough analysis 
of a pancake coil requires taking into account the heat 
transfer between the turns, in addition to that along the 
conductor. Our results for a pancake coil will be pre- 
sented elsewhere. 



FIG. 6: (a) The spatio-temporal evolution of a cryostable 
conductor (« > Kc). The dark (red) spots indicate elevated 
temperatures close to the critical temperature. An individual 
soliton (hot spot) moves by dividing a single temperature peak 
into two twins with subsequent dissipation of one of them. 
The scale on the right and the upper scale correspond to 7 = 
6.7 s"^ and h = 0.4 cm (Eqs. (37) and (38)). (b) The time 
evolution of the temperature profile. Press 'p' to pause/play, 
'q' to rewind, 'f to advance and 'r' to reverse by a frame. 
(The animation works in Windows and Mac Adobe Reader.) 
Alternatively, a snapshot of the temperature profile at jt = 
150. This is the temperature variation along the dotted line 
in Figure |6]^a). 



The practical significance of the dissipative structures 
and solitons forming in the current-carrying wire under 
certain conditions is not immediately apparent. How- 
ever, these are interesting phenomena in their own right 
and it would be worthwhile to try to observe and study 
them experimentally. Their detection can be made by 
conventional electric methods of detection, but a more 
spectacular result could be obtained with the help of the 
real-time magneto-optical imaging [5TJ . 
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